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Possibility of superconductivity from electron repulsion in the Shastry-Sutherland lattice, which 
has a spin gap at half filling, is explored with the repulsive Hubbard model in the fluctuation- 
exchange approximation. We find that, while superconductivity is not favored around the half- 
filling, superconductivity is favored around the quarter-filling. Our results suggest that the Fermi 
surface nesting is more important than the spin dimerization for superconductivity. 
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Superconductivity from electron repulsion has ac- 
quired a renewed momentum from the discovery of high- 
Tc cuprates jj, and the subsequent seminal proposal by 
Anderson Q that the electron correlation should be at 
the heart of the superconductivity. An important, and 
still not fully understood, question dating back to the 
early stage along this line is the relation between the spin 
gap and the superconductivity. A gap in the spin exci- 
tation associated with a quantum spin liquid is, crudely 
speaking, favorable for a singlet-pairing formation. So 
one may naively expect that doping of carriers into a 
spin-gapped Mott insulator will show superconductivity. 
A clear-cut example is the t-J or Hubbard model on the 
two- leg ladder 3], where the un-doped ladder is a spin- 
gapped Mott insulator and the system indeed becomes 
superconducting when doped. 

On the other hand, the spin gap is obviously not a 
necessary condition for superconductivity. The Hubbard 
model on the 2D square lattice or the three-leg ladder 
are examples, which are spin-gapless at half- filling 0, 
but superconduct when dopedQ. 

If we go back to the superconductivity from repul- 
sive electron-electron interactions in a broader context, 
usual understanding is that the effective attraction be- 
tween electrons are mediated by spin fluctuations. @, 
An important difference from the superconductivity from 
attraction is that the effective attraction arises from pair- 
scattering processes across which the BCS gap func- 
tion changes sign, so the effective attraction is wave- 
number dependent and the pairing is anisotropic (typi- 
cally (i-wave as in the cuprates) . Numerically, a quantum 
Monte Carlo calculation that takes care of the relevant 
energy scale 0, || has shown an enhanced pairing corre- 
lation in the repulsive Hubbard model. Analytically, the 
fluctuation-exchange (FLEX) studies which 
is a kind of the renormalized random phase approxima- 
tion based on the Fermi liquid picture, have shown d- 



* Electronic address: |kimuraBcms.phys.s.u-tokyo.ac.jp| 



wave superconductivity with a transition temperature 
T c ~ O(0.01i) for the repulsive Hubbard model. No- 
tably, Tc is two orders of magnitude smaller than t, al- 
though Tc amounts to ~ 100K as in the high-Tc cuprates 
if we take t ~ 0.4 eV. 

Recently, the multiband lattices having disconnected 
Fermi surfaces have been proposed as systems having 
much higher T c 's (> O.Oli) El El El 111 EI While 
the motivation for considering disconnected Fermi sur- 
faces is to raise the Tc, which is lowered down to ~ O.Oli 
in ordinary lattices because the nodes in the BCS gap 
function intersect the Fermi surface: when the Fermi sur- 
face consists of pockets, the nodes can run in between 
the pockets. If we look at the proposed lattices in real 
space, on the other hand, we immediately notice that 
all of them happen to have dimerized structures in some 
way or other. From the spin-gap point of view one might 
consider this reasonable, since the dimerization can favor 
a spin gap when the dimerization is strong. However, we 
cannot identify which of the dimerization and the discon- 
nected Fermi surface contributes to the higher Tc, since 
the dimerization and the disconnected Fermi surface are 
simultaneously satisfied in these lattices. 

Now, there is an important and unresolved question 
of whether the spin gap associated with strong dimeriza- 
tion is a sufficient condition for superconductivity, i.e., 
whether the doped spin-gapped system can always be- 
come superconducting with an appreciable Tc- There are 
several lattices intensively studied from the viewpoint of 
the spin gap in the Mott-insulator phase. Among them 
is the Shastry-Sutherland (SS) lattice (Fig. 1), where 
we have a herringbone array of dimers. This lattice is 
first proposed by Shastry and Sutherland [lj| about two 
decades ago. They found that the ground state in the 
Heisenberg model is spin-gapped when strongly dimer- 
ized, J2/J1 < 0.5, where Ji(Ja) is the spin-spin interac- 
tion within (across) dimer. With the exact diagonaliza- 
tion method for finite clusters, Miyahara and Ueda poj 
have recently obtained a more accurate critical value, 
J%l J\ ~ 0.7, for the appearance of the spin gap. An- 
alytic studies such as a series expansion |2lj . rigorous 
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FIG. 1: The Shastry-Sutherland lattice. 



bounds , and large- iV theories have indicated sim- 
ilar boundaries. 

A recent impetus came from experiments on a cop- 
per compound, SrCu2(B03)2,[24| where the Shastry- 
Sutherland lattice is realized. Experimental results for 
magnetic susceptibility HIS, Cu NQR high-field 
magnetization j24j, ESR l2fl . Raman scatterin g |27l|. in- 
elastic neutron scattering j2|| and specific heat [23] show 
that this compound has a dimerized ground state with a 
spin gap ~ 30K with no long-range magnetic order. [30l| 

So the SS lattice is an appropriate system for the above 
question of whether a spin-gapped system can supercon- 
duct when doped. A recent study [33 for the t-J model 
on the SS lattice has indeed discussed possible supercon- 
ducting transition at low temperatures of O(0. Oli) based 
on RVB type mean field theory j34[ . Usually the SS lat- 
tice is studied with the Heisenberg model, which corre- 
sponds to the half-filled Hubbard model. Here we have 
opted for the Hubbard model, since we can only study a 
finite Hubbard U, but also look at the SS lattice around 
quarter filling, which we propose here to be interesting. 
The band filling controls the shape of the Fermi surface 
after all, so that it should be an important parameter 
for studying the question at hand. Quarter-filled case is 
not unrealistic, since the herringbone structure of the SS 
lattice strongly reminds us of a class of dimerized organic 
crystals, where the band is often quarter filled rather than 
half-filled. 

So here we take the SS lattice to study supercon- 
ductivity in the Hubbard model. We adopt the FLEX 
approximat ion, which has to be extended to four-band 
systems [H [H, H HU to treat the SS lattice that has 
four atoms per unit cell. Superconducting transition has 
been examined with Eliashberg's equation |32j . We shall 
show that the Hubbard model, around half filling, does 
not exhibit superconductivity with significant Tc- So this 
provides an example in which a doped spin-gapped sys- 
tem does not guarantee an appreciable Tc- By contrast, 
the Hubbard model around the quarter-filling exhibits 
superconductivity with a d-wave pairing when the dimer- 
ization is strong t% "C t\. We shall discuss this in terms 
of the shape of the Fermi surface and the pair-scattering 



(a)X (b)|G B | 2 (c)<k 




FIG. 2: Result for \ ( a )> |Gs| 2 (b), and 4>b (c) against k x ,k y 
for the lowest Matsubara frequency for the nearly half-filled 
n = 0.85 with U = 7, ti = 1.8, t 2 = 1.0, and T = 0.025. 



processes on it. 

In the four-band version of the FLEX [3 |H |H EU , 
Green's function G, spin susceptibility x, self-energy £ 
are 4x4 matrices, e.g., Gi m (k,ie p ) with e p being the 
Matsubara frequency. Here I, m refer to the four sites 
in a unit cell, which can be unitary-transformed to band 
indices. We obtain the eigenvalue and the superconduct- 
ing gap function <pi m by solving the linearized Eliash- 
berg's equation |32J. For the susceptibility Xi we quote 
hereafter the value of the largest component when we 
diagonalize its matrix. In the present study, we take 
32 x 32 A;— points and up to 8196 Matsubara frequencies, 
or 64 x 64 A;— points and up to 4096 Matsubara frequen- 
cies. 

We start with the case of near half-filling. Here we 
take the band filling n(=number of electrons/number of 
sites)=0.85, along with the on-site repulsion U — 7 and 
the transfer energy within (across) dimer t\ = ±1.25 
(i2 = 1.0). This value of t\ has been adopted in Ref.|33l] 
as appropriate to SrCu2 (603)2- Corresponding Heisen- 
berg model at half filling has a spin gap because of 
J2/J1 ~ t\lt\ < 0.7. 

The maximum eigenvalue of Eliashberg's equation A 
has turned out to be much smaller than unity (A ~ 0.43) 
at low temperature (0.01 < T < 0.04). We have also 
calculated for t\ = 1.8, which is a favorable value for 
the near-quarter-filled case as we shall see below, but 
A is again very small (A =0.36, 0.37 at T =0.04, 0.01, 
respectively). 

Figure 2 shows the spin susceptibility (x), Green's 
function (G) and the gap function (</>) for the second band 
from the bottom (called B 7 which crosses the Fermi en- 
ergy) for ti = 1.8 and t 2 = 1.0 at T = 0.025. We see the 
susceptibility has no strong peaks, unlike in the square 
lattice which has a large antifcrromagnetic peak in the 
susceptibility around k = (71", 7r) that is relevant to the 
superconductivity. The weak spin structure in the SS 
lattice should be due to the spin gap and/or the spin 
frustration which strongly prevent (s) the long range spin 
correlation when the dimerization is strong|35j. 

If we weaken the dimerization by making t\ sufficiently 
smaller than £ 2 , the system becomes superconductive. 
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quarter filling can be represented by a two-band [37], |40j 
or a single-band 03 |3!J Hubbard model around half 



FIG. 3: The same plot as in Fig. 2 for a weaker dimerization 
with t\ = 0.5. Q represent the nesting vector 
(solid lines). 



Figure 3 shows the spin susceptibility, Green's function 
and the gap function for t\ — 0.5 and t-z = 1.0 at 
T = 0.025, where we obtain a large A = 0.94 close to 
unity (A becomes unity at lower temperature T ~ 0.02). 
In the figure we see the Fermi-surface nesting across 
k x = ±h y is appreciable and the spin susceptibility has 
a strong peak around (0,0), which corresponds to the 
peak around (n, n) on the square lattice (t\ = 0) folded. 
The result indicates that the dimerization (or the spin 
gap) is by no means a sufficient condition for very high 
Tc- If we turn to the very high Tq systems obtained 
in 0, 0, 0, 0, 0|, dimerization causes the discon- 
nected Fermi surfaces accompanied by antiferromagnetic 
spin fluctuations. So the present result is an example in 
which a dimerization works unfavorably for superconduc- 
tivity in a simply-connected Fermi surface. In this sense, 
the disconnected Fermi surfaces rather than the dimer- 
ization is essential for very high T c in [h H li H lul UtI IT^| . 

We have a drastically different situation when we 
change the band filling to quarter filling. Figure 4 shows 
the temperature dependence of A at n = 0.55 with the 
same parameter as those at n — 0.85. We see A is strongly 
enhanced at low temperatures. 

Figure 5 shows %, G and <f> for the lowest two bands 
(called A and B) that cross the Fermi energy for this fill- 
ing. The peak in the spin susceptibility around (0, 0) is 
much stronger than for the half-filling. The spin fluctua- 
tion should mediate the pair-scattering across which the 
gap function has opposite signs, resulting in d-wave su- 
perconductivity. How can this happen when the nesting 
vector is close to (0, 0)? Figure 5(d) depicts the answer: 
the gap function has a d-wave symmetry as in the square 
lattice near half-filling. This is physically natural because 
the SS lattice around quarter filling is effectively a square 
lattice around half filling in the strongly dimerized case 

M»M. 

As mentioned above, it is interesting to compare the 
SS lattice with organic, d-wave superconductors such 
as K-(BEDT-TTF) 2 X.|3g If we assume the dimeriza- 
tion is sufficiently strong, the original system around 
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FIG. 4: For a nearly quarter-filled band (n = 0.55) A is 
plotted as a function of temperature T for U = 7, t\ = 1.8 
and <2 = 1.0. The dotted curve is a least-square fit by a 
fourth-order polynomial. Here we take 32 x 32 fc-points and 
up to 8196 Matsubara frequencies. Error bars for different 
choice of the fc-point mesh and Matsubara frequencies are 
~ 10% at low temperatures. 



filling. |3l| There, superconductivity is enhanced for ~ 
1/4 filled band by strong dimerization as in the present 
study. Although organic materials that can be modeled 
by the SS lattice have not been known, it would be in- 
teresting to search for them. 

In summary, we have studied superconductivity in the 
Hubbard model on the Shastry-Sutherland lattice with 
FLEX. Our analysis shows superconducting transition 
temperature, if any, is very small around half filling de- 
spite the presence of a spin gap due to the dimerization, 
while superconductivity is favored around quarter filling. 
Comparison with the RVB theory for the t-J model |33| 
is interesting because the RVB theory has shown Tc of 
O(0. Oli). While this may naively seem inconsistent with 
our result, the Tc of a conventional RVB theory on the 
SS lattice is much smaller than that (~ O(0. It)) for the 
square lattice[4ll |42| . So the situation about the differ- 
ence between Tc's is similar to this case. 

This work was supported in part by a Grant-in-Aid 
for Scientific Research and Special Coordination Funds 
from the Ministry of Education of Japan. Numerical cal- 
culations were performed at the Supercomputer Center, 
Institute for Solid State Physics, University of Tokyo. 
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FIG. 5: For a nearly quarter-filled band (n = 0.55) x ( a )i 
\G V \ 2 (b), and <f> u (c) (u — A,B) are plotted against k x ,k y 
for the static case (lowest Matsubara frequency) for U = 7, 
ii = 1.8, t2 — 1.0, and T = 0.025. (d) schematically depicts 
the sign of the gap function on the Fermi surface, where Q 
and arrows represent the nesting vector ~ (0, 0). 
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